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Abstract 

The QCD corrections to photon structure functions are defined in a way consistent with the 
factorization scheme invariance. It is shown that the conventional DIS 7 factorization scheme 
does not respect this invariance and is thus deeply flawed. The origins of the divergent 
behavior of photonic coefficient function at large x are analyzed and recipe to remove it is 
suggested. 

1 Introduction 

The recently completed evaluation of order a ^ parton-parton [1, 2] and order aa 2 photon- 
parton splitting functions [3] has confirmed the earlier results [4] based on the evaluation 
of first six even moments. For both the quark distribution functions of the photon and the 
photon structure function, the results (for the latter shown in Fig. 0 exhibit very small 
difference between the next-to-leading (NLO) and next-to-next-to-leading order (NNLO) 
approximations defined in the standard way and evaluated in MS factorization scheme (FS). 

This welcome feature stands in sharp contrast to the large difference, noted already 
in [5], between the standard LO and NLO approximations in the same MS FS. The standard 
formulation of the NLO approximation to photon structure function F^x^Q 2 ) suffers also 
from problems in the large x region where it turns negative. In order to cure these two, 
related, problems, the so called DIS 7 FS has been proposed in [6]. 

The aim of this note is to show that, first, there is no such sharp difference between the 
LO and NLO approximations to F% in the MS FS if these approximations are defined in a 
way consistent with the factorisation procedure. Second, I will argue that the way the DIS 7 
FS is introduced violates the basic requirement of the FS invariance and, consequently, the 
related definition of parton distribution functions of the photon is deeply flawed. To see 
where the problem comes from, I will contrast the definition of DIS 7 FS with the analogous, 
but theoretically well-defined concept of DIS FS in the case of the proton structure function. 

The main source of the difference between the standard treatment of the proton and 
photon structure functions can be traced back to the interpretation of the behaviour of 
parton distribution functions of the photon in perturbative QCD. In [7] I have discussed this 
point at length, but its conclusions have been largely ignored. I will therefore return to this 
point and present additional arguments showing that parton distribution functions of the 
photon behave like a, rather than a/a s as claimed in [3-6] and most, though not all, other 
papers on photon structure functions. I will outline how the theoretically consistent LO, 
NLO and NNLO QCD approximations to F] should be constructed from the results of [3,4], 
Finally, I will discuss the origins of the problems in the large x region and suggest how 
to remove them taking into account their pure QED nature. 


1 



X X 

Figure 1: The comparison of LO, NLO and NNLO approximations to the inhomogeneous 
part of the photon structure function in the MS and DIS 7 FS. Taken from [3]. 

2 Basic facts and notation 


Let us start by briefly recalling the basic facts concerning the various ingredients of pertur¬ 
bative calculations involving (quasi)real photons in the initial state. In QCD the coupling 
of quarks and gluons is characterized by the renormalized colour coupling (“couplant” for 
short) a s (/i), depending on the renormalization scale ji and satisfying the equation 


d a s (/i) 
d In jj? 


PMl*)) 


A) 2/ \ 


ft 

167r 2 


a 3 s (fi) + 


( 1 ) 


where, in QCD with nj massless quark flavours, the first two coefficients, (3q = 11 — 2n//3 
and (3\ = 102 — 38rz//3, are unique, while all the higher order ones are ambiguous. However, 
even for a given r.h.s. of CD there is an infinite number of solutions, differing by the initial 
condition. This so called renormalization scheme (RS) ambiguity 1 can be parameterized in 
a number of ways, for instance by the value of the renormalization scale, usually denoted A, 
for which a s (/i = Ars) = oo. At the NLO the variation of both the renormalization scale 
/a and the renormalization scheme RS={Ars} is legitimate but redundant, ft suffices to fix 
one of them and vary the other. In this paper we shall work in the standard MS RS. 

As nothing in my arguments depends in essential way on the numerical value of f3\, I will 
set in the following (3\ = 0. This assumption, simplifies many otherwise complicated formulae 
and lays bare the essential aspects of the problem. Note, however, that this assumption does 
not amount to working in the leading order of QCD, as all the relevant terms of higher order 
QCD approximations and their relations are kept. 

Factorization scale dependence of PDF of the photon is determined by the system of 
inhomogeneous evolution equations 


d E(x,M) 

d In M 2 
dG(a;, M) 

d In M 2 
dq NS (x,M) 
din M 2 


+ Pqq <8> S + 2 UfP qG <g> G, 

(2) 

k G + P Gq ® E + P GG ® G, 

(3) 

&NS + -PnS <8> <?NS, 

(4) 


1 In higher orders this ambiguity includes also the arbitrariness of the coefficients R, i > 2. 
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where 


n f n f 

E(x,M) = ^2qf(x,M) = '^2[q i (x,M) + q i (x,M)], 

i =1 i= 1 

n f 

q m (x, M) = (e- - (e 2 )) {q^x, M ) + ^(x, M)) , 

2=1 

&ns = ^ns = 6n/ ((e 4 ) — (e 2 ) 2 ) , fcs = fakg] 5y, = 6n/(e 2 ). 

The splitting functions P tl and k t are given as power expansions in a s (M): 


ki(x,M) 


Pij(x,M) 


a 

2tt 


kf\x) + 


a s (M) 
2t r 


k^\x) + 


cts{M) 

2tt 


kf\x) + • • • 


q a (M) 

2n 




(5) 

( 6 ) 

( 7 ) 

( 8 ) 

(9) 


where i = g, G, NS and ij = qq, qG, Gq, GG, NS. The leading order splitting functions 
kf\x) and P^\x) are unique , while all higher order ones k q \ kg , Pjfi\j > 1 depend on the 
choice of the factorization scheme (FS). According to the factorization theorem, the photon 
structure function Ff(x, Q 2 ) is given as the sum of the convolutions 


1 

x 


Ff {x, Q 2 ) 


q NS (M) 0 C q (Q/M ) + S NS C^ + (10) 

(e 2 )E(M) 0 C q {Q/M ) + (e 2 )S E C 7 + (e 2 )G(M) 0 C G (Q/M) (11) 


of PDF and coefficient functions C q , C G , C 7 admitting perturbative expansions 


C q (x,Q/M) 
C G (x, Q/M) 
C 7 (x, Q/M) 


5(1-x) + Q/M) + ■ • •, 

277 


where the standard formula for C\° 


reads 


Cf\x,Q/M) = (x 2 + (1 - x ) 2 ) ln AID—1 + 8xfl -x) - 1. 


( 12 ) 

(13) 

(14) 

(15) 


1 am using the terms “renormalization” and “factorization” scales in standard way as dis¬ 
cussed, for instance, in [8-10]. However, some authors use these concepts in a very different 
way. For instance, in [11] the renormalization scale denotes the argument of the couplant 
a s in the expansion © of the splitting functions, rather than in the expansions mm ° f 
the coefficient functions as in [8-10] and the present paper. This alternative definition of 
the renormalization scale, though mathematically legitimate, lacks physical motivation as it 
fails to make the crucial difference, emphasized long time ago by Politzer [12], between the 
ambiguities of the treatment of perturbatively calculable short distance physics, embodied 
in the renormalization scale and described by the coefficient functions, and those of the large 
distance, nonperturbative effects that go into the PDF and induce their dependence on M. 
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Figure 2: Diagrams defining the pointlike part of non-singlet quark distribution function of 
the photon in the leading logarithmic approximation. 

The renormalization scale /i, used as argument of aq(/i) in (H2H14j) is in principle inde¬ 
pendent of the factorization scale M. Note that despite the presence of a s (p) in (fT2HT4l) . the 
coefficient functions C q , Cq and C 1 are, if calculated to all orders in a s , independent of p as 
well as of the RS. On the other hand, PDF and the coefficient functions C q , Cq and C 7 do 
depend on both the factorization scale M and factorization scheme, but in such a correlated 
manner that physical quantities, like P7, are independent of both M and the FS, provided 
expansions flEHHD and (TH2HTT1) are taken to all orders. In practical calculations based on 
truncated forms of mm and mm this invariance is, however, lost and the choice of both 
M and FS makes numerical difference even for physical quantities. 




3 Parton distribution functions of the photon 


The general solution of the evolution equations (US) can be written as the sum of a particular 
solution of the full inhomogeneous equation and the general solution of the corresponding 
homogeneous one, called hadronic. In rest of this note 1 will for technical reasons restrict the 
discussion to the nonsinglet quark distribution function of the photon and the corresponding 
nonsinglet part m of photon structure function. To simplify the formulae 1 will also drop 
the subscript “NS” and set <5 ns = 1 everywhere. 

As is well known, the subset of the solutions of the evolution equation 0]) with the 
splitting functions including the first terms k ^ and P(°) only and vanishing at some scale 
Mo results from the resummation of the contributions of the diagrams in Fig. ((21). These, 
so called pointlike solutions, which start with the purely QED vertex 7 —> qq , define the 
standard “leading order” approximation and have, in momentum space, the form (again 1 
will drop the specification “pointlike” throughout the rest of this paper) 


q(n, M 0 , M) 


47T 

a a (M) 


a{n ) 


1 


a,(M) y- 2 pW W/ft 


(16) 


where 

/ \ _ a k^°\n) 

“ (n) = 2 tt/3 0 1 - 2P<»>(n)//V ' ) 

As argued in [7] the fact that a s (M ) appears in the denominator of (fTfil) does in no way 
mean that q oc a/a s . It is obvious that provided M and Mo are kept fixed when QCD is 
switched off by sending A —> 0 the expression m approaches 

q(x, M, M 0 ) —* ^£; ( 0 ) ( x) In ^ > ( 18 ) 

corresponding to purely QED splitting 7 —> qq. Note that this limit holds even if we include 
in 0 the term proportional to k^\ which gives the lowest order QCD contribution to the 
inhomogeneous splitting function k q . 
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To see what is wrong with the conventional way of arriving at the claim that q oc a/a s 
let us recast the evolution equation 0 into the equivalent form 


d q(n,Q) _ 47 T 

da. 


a k^\n) a P^(n) . _ 

+ x- —q(n,Q) + 


(3q \_2n a 2 ' 2n a s 
and take into account just the first term on its the right. Trivial integration then yields 

a 47 t Ad°)(n' 


(19) 


q(n,Q) = 


2 tt f3 0 at a {Q) 


A, 


( 20 ) 


where A denotes arbitrary integration constant specifying the boundary condition on the 
solution of m- Choosing A = 0 might, but should not, mislead us to the usual claim that 
q oc l/a s , because (ITU is equivalent to 0 , which in our approximation of keeping just the 
first term on its r.h.s. contains no trace of QCD, being of purely QED nature! Taking the 
difference of (EDI) for Qi and Q 2 and inserting the explicit expression for a s (M) we arrive at 


q(n,Qi) - q(n,Q 2 ) = f 


oisiQl) a s (Ql 


= -^-k (0 \n) ln^i, 
V 2 


( 21 ) 


the purely QED expression we can get directly from 0 . We can further recast m into the 
form similar to (USD 


q{n, Q 1 ) - q{n, Q 2 ) = 


47r 


a 


k^\n) 


a s (Ql) 2ir/3 0 d 


1 - 


(XsiQV, 

<*s(QZ 


( 22 ) 


with d — 1. Including also the lowest order QCD term proportional to P^°\n) modifies the 
result slightly, but non-essentially by replacing d = 1 with d — 1 — 2P^°\n) / (3q. Note that 
since Pjfq (0) = 0 the integral of quark distribution function is actually unchanged by QCD 
effects, but the standard claim would still be that it “behaves” as 1 /a s . 

By the same reasoning we could “prove” that, for instance, also the vacuum polarization 
n(Q 2 ) “behaves like l/a s ”! Indeed, taking the derivative of II(Q 2 ) with respect to InQ 


D{Q ) = 


2\ _ dn(Q 2 ^ 


dn(Q 2 )da s (Q) 


dlnQ 2 da s (Q) dlnQ 2 
which defines the Adler function, and rewriting it in the way analogous to m we get 

dn(Q) = D(Q 2 ) = 47T l + Er=i d k a k s (Q) 

d a s (Q) da s (Q)/d In Q 2 

Trivial integration then yields the advertised conclusion 

47T 1 


(23) 


a 2 + 


n(Q 2 ) = 


(3o ol s (Q) 


+ 


(24) 


(25) 


Let me finally add another, cpiite different, argument demonstrating that PDF of the photon 
behave as O(a). Consider the Mellin moments of the nonsinglet part of F 7 (x) = F^ s (x, Q)/x 


F J (n, Q) = q(n, M)C q (n, Q/M) + C 7 (n, Q/M ) 


(26) 
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for the general pointlike solution of ED, which like (1TT)1) vanishes at M = M 0 . As F 7 (n, Q) 
is independent of the factorization scale M, we can take any M to evaluate it, for instance 
just M 0 . However, for M = Mq the first term in (12611 vanishes and we get 2 


a 


F ^Q) = s 


C^HQ/Mo) + ^-C^iQ/Mo) + QA>) Cf(Q/M„,Q//j) + 


(27) 


i.e. manifestly the expansion in powers of a s (/i ) which starts with O(a) pure QED contri¬ 
bution (a/27r)C^°' > and includes standard QCD corrections of orders a k s , k > 1. Clearly this 
expansions vanishes when QCD is switched off and there is no trace of the supposed “a/a s ” 
behaviour. I will come back to this expression in Section 4. 


4 Defining LO and NLO approximations for F 7 (x, Q 2 ) 


Although semantics is a matter of convention, I think it is wise to define the terms “leading”, 
“next-to-leading” and higher orders in a way which guarantees that they have the same 
meaning in different processes. Recall that for the case of the familiar ratio 


R e + e -(Q) 


cr(e + e —*• hadrons) 
cr(e + e _ —> /U + /i _ ) 



(1 +r(Q)) 


(28) 


the prefactor ‘2>Y^i=i e ‘ii which comes from pure QED, is usually subtracted and only the 
QCD correction r(Q) is considered for further analysis. For the quantity (1281) the terms 
“leading” and “next-to-leading” thus apply only to genuine QCD effects as described by 
r(Q), which starts as a s / tt. Unfortunately, this practice is ignored in most analyses of F 7 , 
which count the QED term as the “leading order” [3-6]. 

In [7] I have proposed the definition of QCD approximations of F 7 which follows closely 
the convention used in QCD analysis of quantities like EHP. It starts with writing the (point¬ 
like nonsinglet) quark distribution function q(M) as the sum of the purely QED contribution 


qqed(M) = In YJ 2 


2tt 


Mn 


(29) 


and the QCD correction satisfying the inhomogeneous evolution equation 

d<?QCD (M) 


din M 2 


2n 


^-fc (1) + P (0) q qed( m) 

Z7T 


a. 




2 r 


£fc( 2 )+P« gQED ( M) 


+ 


+ ^e (0) 9QC D (M) + (g) P (1 > 9 qcd(M) + 


(30) 


The latter differs from that satisfied by the full quark distribution function not only by the 
absence of the term (a/2ir)k ( '°' > but also by shifted appearance of higher order coefficients 
A 77 ; i > 1. For instance, the inhomogeneous splitting function k\ ^ enters (1601) at the same 
order as homogeneous splitting function P and thus these splitting functions will appear 
at the same order also in its solutions. Similarly, the simultaneous presence of kq 2 ' 1 and P (l> 
in the 0(a 2 ) term of the inhomogeneous part of (1301) implies that the NLO QCD analysis 

2 In the rest of this paper the dependence on the Mellin moment variable n will be suppressed. 
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of F 7 requires the knowledge of k^ etc. In Table |T| the terms included in the standard 
definition of the LO and NLO approximations are compared with those corresponding to my 
definition of these approximations. The difference between the two definitions is substantial. 



standard definition 

my dehnition 

QED 

does not introduce 


LO QCD 

Jfe(°), P(°) 

c!f\ fed), C\ 1] , P (0) 

NLO QCD 

fc(°), P (0) , fc (1) , C^\ 

fcW, C^\ p(°), k {2 \ c!f\ P« 


Table 1: Contributions included in the standard and correct definition of LO and NLO. 


5 Factorization schemes and their choice - the proton 


Let us first recall the definition of factorization schemes in the case of the (nonsinglet) proton 
structure function. Denoting F p = F^ s /x we can write its moments as the product 

F p (Q) = q(M)C q (Q/M) (31) 

of the coefficient functions C q (Q/M), given in (11211 . and the quark distribution function q(M) 
of the proton. If both the homogeneous splitting function © and coefficient function m 
are calculated to all orders in a s , the product on the r.h.s. of en is independent of both 
the factorization scale M and the factorization scheme. It must thus hold 


d F(Q) 

d In M 2 


= q(M) P(M)C q (Q/M) + C g (Q/M) 


d F{Q) 


d a 


U) 


= 0 


(32) 


If, however, these expansions are truncated, the resulting finite order approximations for 
F(Q) will depend on both the factorization scale and scheme. For a finite order approxima¬ 
tion to be theoretically consistent, its dependence on all free parameters must be formally 
of higher order in a s than those included in the approximation. From (|S2D we get, denoting 
f(M) = d/(M)/dlnM 2 , 

P(M)C q (Q/M) + C q {Q/M) = 0 (33) 

and expanding the splitting and coefficient functions to order a s we arrive at the relation 
C^\Q/M) = -P (0) => CW(Q/M) = P {0) In(Q 2 /M 2 ) + C^l). (34) 

Taking into account the first two terms in © the solution of © reads 

_2p(°) 

2 Pd) a s (M) 


q{M)=A(a s {M)) ^ 


exp 


(35) 


A) 27r 

where the (n-dependent) constant A specifies the boundary condition 3 . Inserting (HTHl) into 

(EH) yields 

_ 2p(°) 

2 Pd) a s (M) 


F P (Q) = A (a s (M)) A) 


exp 


A) 


2vr 


1 + ’(<?/") 


(36) 


3 We can also use the above expression for M and Mq and by forming their ratio obtain instead of 03 
the usual relation between q(M) and q(Mo). 
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and expanding the exponential we get 


-2p(°) 


(37) 


F-(Q) = 31 (a.(M)) A> 1 + ^ ^ C W(Q/M) - 
FS invariance of (J36l) implies that the non-universal functions Cq 1 * and Pq\ } are related 

Cf(l) = ^ + s «■ P (1 > = § (Cf(l) - K ) , (38) 


where the quantity k = /t(n) is factorization scheme invariant . In other words either P ( d 
or Cq \ but not both independently, can be chosen at will to specify the FS. For instance, 
F P (Q ) can be written as a function of Cq* explicitly as 


_2P(°) 

F P (Q) = A(a s (M)) A) exp 


(c ,i, _ k) 




(39) 


The above expression is independent of Cq 1 ' ) to the order considered but not exactly and 
thus its numerical value does depend on the choice of Cq 1 * . Two points are worth noting. 

First, the cancelation mechanism based on the relation (1381) operates independently of the 
value of a s as well as for any fixed boundary condition specified by the constant A = A(n). 
These constants, being determined by the behaviour of q(M ) at asymptotic values of M, 
provide unambiguous way of specifying the initial condition on the solution of the evolution 
equations. The same boundary condition must therefore be used in all FS. If the boundary 
condition on q(M ) is specified by the value q(M 0 ) at some initial M 0 , the situation is slightly 
more complicated. Imagine, we have q(M, M 0 ) given as the solution of Q in a FS specified 
by P^\ or equivalently Cq 1 * , and with the boundary condition given by q(M 0 ) at some M 0 . 
To get the boundary condition in a FS specified by P^\ or equivalently C^\ which would 
yield the same proton structure function m we first use (TIH to convert the information on 
g(Mo) into the knowledge of the constants A = A{n). Once we have them we can again use 
(ESI) to compute the appropriate boundary condition q(M 0 ) in the new FS: 


q{M 0 ) = q(M 0 ) exp 


2 pb) - p 




a s (M 0 


A) 


2t r 


= q(M 0 ) exp 


(c (1) - c ,(1) ) 


—(i)\ a s (M 0 ) 


2 tt 


(40) 


This boundary condition is different from the one in the original FS at the same Mo, and 
so would also be the solution q(M, M 0 ), but after inserting into (HTTl) we get the same F P (Q ) 
since the change of q(M,M 0 ) is compensated by the accompanying change of C q (Q/M). If 
we perform an analysis of experimental data on F P {Q ) by fitting the parameters specifying 
the initial q(M 0 ) (and Aqcd) we can thus work in any FS and choose in principle any initial 
M 0 and get the same results 4 for F P (Q). 

Second, the expression m connects quantities, Cq 1 '* and Pd), which come from evalua¬ 
tion of Feynman diagrams at different fixed orders: the former at 0(a s ) order whereas the 

4 The independence of the chosen FS, as well as of the value of Mo does, however, requires sufficiently 
flexible form of the parametrization of q(Mo). 





















latter at 0(0%). This is due to the fact that the leading-logarithmic terms proportional to 
powers of P d) are resummed to all orders thereby generating the first term on the r.h.s. of 
©. Only after this resummation does Cq 1 * contribute to F p at the same order as Pd). 

In the standard MS FS both Cq^(l) and P (V) are nonzero. From the infinity of other 
possible schemes only the so called “DIS” FS is regularly used. In this scheme one sets 
Q = M and Cq ^ = 0 in order to keep at the NLO the same relation between the proton 
structure function F(Q ) and the quark distribution function q(Q ) as at the LO. In some 
sense opposite to the DIS FS is the FS Pd) — o. In this FS the evolution equation for quark 
distribution function has the same form at the NLO as at the LO. As noted in [13] this FS 
corresponds 5 to the point of local stability [14] of the expression considered as function 
of M, ji and Cq 1 ' 1 . 


6 Factorization schemes and their choice - the photon 

For the hadronic part of the quark distribution function of the photon the factorization 
operates in exactly the same way as for the proton. Consequently Cq 1 '* can again be used to 
label different factorization schemes and also the relations do hold. 

For the pointlike part of the quark distribution function of the photon the situation 
is slightly more complicated as the expression for photon structure function P 7 = F^ s /x 
involves another coefficient function, namely the photonic C 1 

F\Q) = q(M)C q (Q/M) + C 7 (Q/M). (41) 


Consequently, the factorization scale invariance implies 


F\Q) = q(M)Cq(Q/M) + q(M)Cq(Q/M ) + C^Q/M) = 0 

= P(M)C q (Q/M) + C q (Q/M)} q(M) + k(M)C q (Q/M) + C^Q/M) (42) 


The expression in the square bracket vanishes for the same reasons as for the hadronic 
structure function and we are thus left with 


P 7 (Q) = 


a 

2tt 


<! fc (0) + C^*(Q/M) + ^ C {1 *(Q/M) + fed) + fenc'd) (Q/M) 


2tt 


Jfc(°) + ^Mfcd) + ■ ■ ■ 

2n 

C<®(Q/M) + ^-C^(Q/M) + 

CXc 


! + ^Cf(Q/M) + ... 


+ 


(43) 


2tt 


+ 


= 0 . 


The argument of a s in the expansion of the coefficient functions C q and C 7 , i.e. the renor¬ 
malization scale /i, is in general different from the factorization scale M. However, I did not 
write out the argument of a s in the last line of © because the coefficient standing by it is 
independent of it. The Eq. © implies 

CW(Q/M) = -k i0 \ (44) 

C^iQ/M) = -fed) -fenc'd)(Q/M), (45) 

5 This holds exactly for /3i = 0, for realistic values of /?i the saddle point is slightly away from pd) = 0. 
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which upon integration and taking into account the relations (Id411 and (El) yields 


C^iQ/M) 

C?(Q/M) 


fc (0) hi(Q 2 /M 2 ) +Gf (1) 


fc(°) p(°) 
2 


In 2 (Q 2 /Af 2 ) + (fc (1) + 4 (0) cy (1)) ln(Q 2 /M 2 ) + Cpl). 


(46) 

(47) 


Let us now consider the pointlike solution which results from taking into account beside the 
pure QED splitting function k ® the lowest order QCD splitting functions P^ and k^: 


q(M,M 0 ) = 


47ra 

(c(M)V- 2pm C 

a 

j (a,(M)Y 2Pm C 

Ad 1 ) 

a s (M) 

\a s (M 0 )J 

2 t r 

\a s (M 0 ) J 

p( o)’ 


(48) 


where a = a{n) is defined in (El). This solution satisfies q(Mo, Mo) = 0 and for small a s 
behaves as 


a 


q(M,M 0 ) — —L m 
Ztt 




2tt V 2 

Multiplying © with C q and adding the lowest two terms of C 1 we get 


T M2 
• Lm= mJ' 


(49) 


2vr 


—F\Q) = L M k (0) + Cf (Q/M) 


a 


(50) 


+ 


27T 


C^(Q/M) + (fc^ + k( 0) CP(Q/M)) L m + 


&(°) p(o) 


r 2 

Ll M 


Replacing Cy 1 with El) and C q ' with analogous expression (I38[) . we hnd that all terms in 
the square bracket dependent on M cancel out 6 and we are left with expression which is a 
function of Q/M 0 only. Moreover, as expected from Ell . it is just C^\Q/Mq)\ 

Note that the only ambiguity in m comes from the freedom in the choice of the ex¬ 
pansion parameter a s (/x), i.e. the choice of the renormalization scale /j and renormalization 
scheme. Because these ambiguities appear first at order ex 2 , the lowest order QCD correction 
given by C^\Q / M 0 ) must be unique! Consequently, C? ?1) as well as the combination 

fc (1) + fc (0) C'W(l) = k (51) 

must also be unique, i.e. factorization scheme independent. The above relation implies that 
the Erst non-universal inhomogeneous splitting function Ad 1 ) is, similarly as Pd) ; a function 
of C^(l) (or the other way around). So Cq 1 \ 1) and its higher order partners Cq } can be 
chosen to label the factorization scale ambiguity also for P 7 . Once these quantities (which 
are functions of n, or equivalantly, x) are chosen all other free parameters are determined 
by relation like those in (El) or ED. 


7 What is wrong with the DIS 7 factorization scheme? 

The freedom in the definition of quark distribution functions of the photon related to the 
non-universality of Cq\j > 1 is intimately connected with higher order QCD corrections. 
The purely QED contribution, given by the sum of first two terms in is analogous 

6 The analogous cancelation in the pure QED contribution to F 1 is trivial. 
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to QED contribution 3 Y^i=i e i i n and represents an input into QCD analysis. It is 
manifestly of order a° = 1 , but since in the standard approach the quark distribution 
function is - incorrectly as I have argued - assigned the order 1 /a s , the lowest order QCD 
contribution appears to be of lower order than . The latter is consequently assigned 
to the “next-to-leading” order [3, 5, 6 ] and treated in a similar way as the lowest order 
QCD coefficient function C^. This procedure has been motivated in part by the fact that 
C-y°\x, M) as given in (1151) turns negative at large x and even diverges to — oo when x —> 1. 
This has lead the authors of [5] to introduce “A factorization scheme avoiding the common 
perturbative instability problems encountered in the large-x region. ’’ In this so called DIS 7 
factorization scheme "... we can use the same boundary conditions for the pointlike LO and 
HO distributions 

Opl( x iQo ) = 1pl( x iQo) = Gf L ( x ,Q 0 ) = 0 ( 52 ) 

without violating the usual positivity requirements. ” [5]. 

Their procedure amounts to of redefining the quark distribution function of the photon 
by absorbing the pure QED term C^ in q(M,M Q ) according to (see eq. (5) of [ 6 ]) 

q(M, M 0 ) = q(M,M 0 ) + |-Cf (1). (53) 


As this redefinition involves the QED quantity Cj , and the notion of quark distribution 
function inside the photon is well-defined also in pure QED, the above procedure must make 
sense even there, without any QCD effects. Moreover, as the QCD contribution depends on 
the numerical value of a s it cannot cure any problem of the pure QED part. 

However, in QED it is straightforward to see that in getting rid off the troubling C' 7 °' ) term 
the procedure proposed in [5] violates the requirement of factorization scheme invariance. 
Recall that in pure QED the contribution to F 7 coming from the box diagram regularized 
by explicit quark mass m q reads [15] 


foEDW) = Ac { °\Q/m.) = 


QED 


a 


2 vW 7 v ~ w " q ' 2 vr 


fc (0) ln-^ + Cf (1) 

mi ' 


Introducing the arbitrary scale M, we can split it into 

<?qed (M) = hi 

2tt m z q 


(54) 


(55) 


interpreted as quark distribution function, and CP(Q/M) given in IIIQlt : 

*3ed«) = 9qed(M) + |-iGf(Q/M ). (56) 

The sum (ES|) is manifestly M- independent. We can now redefine gQED(Af) by adding an 
arbitrary function / = f(n) according to 


QL 

QQED,f(M) = qQEp>(M) + —f. 

In order to keep the sum 

^edW) = 9QED,/(M) + f-C^Q/M) 


(57) 


(58) 
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independent of M and / and equal to 0 induces the correlated change of C- 


(0) 


C ( °}{Q/M) = Cf°\Q/M)-f. (59) 

Note that the DIS, factorization scheme of [5] corresponds to / — cf (1)- We can write 
down the evolution equation for gq ED j(n, M) in the “/-factorization scheme” 


d<?Q E p,/(M) _ _a_, (o) 

d In M' 2 2ti ’ 


(60) 


which is the same for all gq ED j(M). What we are, however, not allowed to do is to use the 
same boundary condition for all g QED j(M), i.e. precisely what has been assumed in (1521) ! If 
we do that and impose the boundary condition g QED j(M = m q ) = 0, we get 

= ( 61 ) 

which depends on the choice of / and only for / = 0 coincides with the correct result 0 . 
To get the latter we must impose on the solution of the evolution equation (EDI) in the /- 
factorization scheme the appropriate boundary condition: gq E o,/(M = m q ) = (a/27r)/ . The 
reason the QED expression (P|) depends on the FS specified by the function / = f{n) is 
clear: there is no all order resummation of the LL and NLL terms describing the multiple 
photon emission off the qq that would give rise to terms analogous to the first and second 
terms on the r.li.s. of ( 0 - 

In QCD the redefinition (1571) (with gq EE replaced with full quark distribution function 
q(M,M 0 )) induces not only the shift (j551) of C^ but also those of k O and 


kf ] = k {1) - P (0) / (62) 

C%(Q/M) ee CV{Q/M)-C£\Q/M)f. (63) 

As in pure QED the DIS 7 factorization scheme corresponds to / = C 7 °\l). In [5] and most 
other NLO QCD analysis the term proportional to C 7 1 ' ) is assigned to the NNLO and thus 
the above second relation 0 is not written out explicitly. 

The above substitutions (I59H63I) are legitimate but as in QED what we are not allowed 
to do is to impose on q/(M, Mq) the same boundary condition for all /. If we do that we 
straightforwardly find that the photon structure function can be written as 


F(Q) 


a 

2tt 


(Cf(Q/M 0 ) - /) + ^ (C«(Q/M„) - Cf (Q/M)/) + • • ■ 


(64) 


which coincides with 0 only for / = 0. Even if we sum 0 to all orders of a s the result 
does depend on the chosen /, violating the fundamental requirement of factorization scheme 
independence of physical quantities!! As already emphasized, the QCD corrections cannot 
compensate the /-dependence of the first pure QED term. In fact even the all order sum of 
QCD terms in 0 is /-dependent. 


8 Removing the problem in large x region 

Because the function C^°\x,Q/M) is of pure QED origin, one might expect it to cause 
problems also for leptonic structure function of the photon. This structure function, which 
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had been measured at LEP, is given by the same expression (USD as for quarks, except for 
the replacement of m q with the lepton (muon or electron) mass mi: 

^r F 7ept( x iQ) = + (! - X ) 2 ) In Q ^ X - + 8 x(l -x)-l 

Lx. I f Jb 

W 2 

= (x 2 + (1 — x ) 2 ) In -j—jj- + ( x 2 + (1 — x ) 2 ) In 4 + 8 x(l — x) — 1. (65) 

The lepton mass regularizes mass singularities of (©), as does the quark mass in (©)• Since 
in the kinematically accessible region W 2 > 4 mf, corresponding to x < 1/(1 + 4 m 2 /Q 2 ), 
both the function B(x) = ( x 2 + (1 — x) 2 ) In4 + 8 x(l — x) — 1 and the first term in (1651) 
are positive, there is thus no problem with positivity of F^ ept . The same holds for the QED 
contribution to quark distribution function with m q as infrared regulator. 

To identify the kinematical regions contributing to various parts of © or (© we return 
to the x-space, where the limits on the virtuality r (see Fig. n of the quark included in the 
definition of the pointlike quark distribution function of the photon are given as [17] 

m l . . Q 2 _ 

^rnin ^ T ^ — T'max* (OOJ 

1 — X X 

The potentially troubling term ln(l — x) in (fh5l) comes from the lower limit on r. So long 
as it appears in the combination ln(Q 2 (l — x)/m 2 x) as it does in ([1711) it is harmless as 
the whole term ln((5 2 (l — x)/m q x) stays, as argued above, positive. The problem, does, 
however, arise when we replace m q in this combination with an arbitrary initial M 0 , which, 
moreover, is typically around 1 GeV and consider (HuEl) for x > 1/(1 + 4Mq/Q 2 ), where 
ln(Q 2 (l — x)/Mqx) < 0. In this case we should drop the term proportional to ln(l — x) 
which comes from the region cut-off by the introduction of Mq. Moreover, we should also 
drop all other terms which come predominantly from the lower integration range in m- As 
argued in [7,17] this implies that Cy 1 as given in (1151) should be replaced with 

C^ 0 ) (x, Q/M) = \x 2 + (1 - x) 2 ] In + \x 2 + (1 - x) 2 ] In - + 6a;(1 - x) - 1, (67) 

M z x 

which implies decent behaviour of 

(a;, 1 ) = [a ; 2 + (1 — x) 2 ] In — + 6 x(l — x) — 1 ( 68 ) 

r l J x 

when x —> 1. Note that the expression (1HH1) coincides with that used in [18]. 

9 Summary and Conclusions 

The standard definition of LO and NLO approximations to photon structure function are 
shown to be inconsistent with the requirement of factorization scale and scheme invariance. 
The origin of this shortcoming is traced back to the usual but incorrect claim that quark 
distribution functions of the photon behave as a/a s . 

Theoretically consistent definition of LO, NLO and NNLO approximations for the photon 
structure function are constructed and the potentially troubling behaviour of C\°\x,Q/M) 


13 





at large x is removed by discarding the terms coming from the region outside the validity of 
perturbation theory. 

The so called DIS 7 factorization scheme invented in order to cure the mentioned problem 
with large-x behaviour of C^\x,Q/M) is shown to be ill-defined as it is based on the pro¬ 
cedure which violates the factorization scheme invariance of the photon structure function. 
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